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Using Wilson renormalization group, we show that if no integrated vector operator of scaling 
dimension —1 exists, then scale invariance implies conformal invariance. By using the Lebowitz 
inequalities, we prove that this necessary condition is fulfilled in all dimensions for the Ising uni¬ 
versality class. This shows, in particular, that scale invariance implies conformal invariance for the 
three-dimensional Ising model. 


I. INTRODUCTION 

Conformal symmetry plays a considerable role both in 
high energy and condensed matter physics. There has 
been a renewed interest in recent years, particularly be¬ 
cause of the AdS/CFT conjecture [lj and the successful 
use of conformal methods in three-dimensional critical 
physics The groundbreaking papers of the 1970s 

and 1980s l8l 1 13l| solved two fundamental issues in two 
dimensions: First, scale invariance implies conformal in¬ 
variance under mild assumptions [12i . Il3j | and, second, 
conformal symmetry enables us to solve most of the scale 
invariant problems that is, to determine critical expo¬ 
nents and correlation functions Q 

An important ingredient for the exact solution of two- 
dimensional conformal models is the existence of an infi¬ 
nite number of generators of the conformal group. In 
higher dimensions, the number of generators is finite, 
and we could naively conclude that symmetry arguments 
alone are not sufficient to solve a model in the critical 
regime. However, it is well-known that scale-invariant 
theories are in a one-to-one correspondence with the fixed 
points of the Wilson Renormalization Group (RG) ITH . 
and that the fixed point of a theory completely deter¬ 
mines all the correlation functions of a critical model for 
sufficiently small wavenumbers. Therefore, at the level 
of principles, scale (and a fortiori conformal) invariance 
is sufficient to determine all the universal critical prop¬ 
erties of a model. Of course, in practice, the computa¬ 
tion of these critical properties requires us to solve the 
functional Wilson RG equations. This is a formidable 
task that we do not know how to carry out without ap¬ 
proximations. Any supplementary information, even if 
redundant, is therefore welcome and this is what con¬ 
formal invariance could provide. A breakthrough in this 
direction has been achieved these last years with the con¬ 
formal bootstrap program @-0|, which led to the exact 
(although numerical) computation of the critical expo¬ 
nents of the Ising model in three dimensions assuming, 
among other things, conformal invariance. 

In parallel, a large activity has been devoted to under¬ 
standing the relation between scale and conformal invari¬ 
ance in - or close to - four dimensions. It has been proven 
to all orders of perturbation theory that scale invariance 
implies conformal invariance fid ] in four-dimensional uni¬ 


tary and Poincare invariant theories. Moreover, there are 
strong indications that a non-perturbative proof could be 
at reach in this dimension 111 71-1 19f. 

Despite decades of effort, it is still an open question to 
know whether a typical statistical model is conformally 
invariant at criticality in three dimensions. The aim of 
this article is twofold. First, we derive a sufficient condi¬ 
tion which, when fulfilled, ensures that scale invariance 
implies conformal invariance. In the second part of the 
paper, we prove that this condition is fulfilled in any di¬ 
mension for the euclidean Li model. 

The rest of the paper is organized as follows. In 
Sect. El we make a brief review of the nonperturbative 
renormalization-group. We then recall in Sect. 1III1 the 
relation between scale invariance and the NPRG. By us¬ 
ing the same methods, we generalize these considerations 
to the case of conformal invariance in Sect. ED In Sect. 
El we finally derive a sufficient condition for the valid¬ 
ity of conformal invariance in scale invariant models. We 
show, on general arguments that this condition is ex¬ 
pected to be fulfilled in O(N) models (and in generaliza¬ 
tions thereof). In Sect EH these considerations are made 
rigourous for the Ising universality class. We give our 
conclusions in Sect IVIII 


II. NONPERTURBATIVE 
RENORMALIZATION-GROUP FORMALISM 

The proof of conformal invariance in all dimensions 
presented below is intimately related to the deep struc¬ 
ture of the Wilson RG 1 and scale invariance. We there¬ 
fore start by recalling, in the case of the Z 2 model, 
the formalism of the modern formulations (sometimes 
called Nonperturbative RG, or functional RG) of the Wil¬ 
son RG [22M26I] . The coarse-graining procedure at some 
RG scale k is implemented by smoothly decoupling the 
long-wavelength modes p[\q\ < k) of the system, also 
called the slow modes, by giving them a large mass, 
while keeping unchanged the dynamics of the short- 


1 The history of the relation between conformal invariance and 
Wilson RG is long, see for example f2Cj|, DTI . 
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wavelength/rapid ones ip(\q\ > k). This decoupling is 
conveniently implemented by modifying the action or the 
Hamiltonian of the model: 5[(^] —» S[y>] + A5fc[y>] where 
AS//^] is quadratic in the field and reads, in Fourier 
space, ASfcjy] = 1/2 f q R k (q 2 )ip(q)ip(-q). The precise 
shape of R k (q 2 ) does not matter for what follows as long 
as it can be written as 

Rk(q 2 ) = Z k k 2 r(q 2 /k 2 ) (1) 

where Z k is the field renormalization factor and r is a 
function that, (i) falls off rapidly to 0 for q 2 > k 2 - the 
rapid modes q>{\q\ > k) are not affected by A S k - and (ii) 
goes to a constant for q 2 = 0 - the slow modes <p{\q\ < k) 
acquire a mass of order k and thus smoothly decouple. 
The partition function, which now depends on the RG 
parameter k, reads: 

Z k [J} = J Vip exp ^-S[q>] - AS k [ip] + J Jtpj (2) 

where the field J is a source term which corresponds to 
the magnetic field in the Ising model and where the ultra¬ 
violet (UV) regime of the functional integral is assumed 
to be regularized at a momentum scale A, see for instance 
[2?t for a lattice regularization in this formalism. It is 
convenient to define the free energy W k [J] = In Z k [J] 
and its (slightly modified) Legendre transform by 

r fe [</] + W fc [J] = f j R k {\x-y\)<j>{x)<j>{y), (3) 

with <f>(x) = SW k /SJ(x), R k (\x — y\) the inverse Fourier 
transform of R k {q 2 ) and where the last term has been 
added for the following reason. When k is close to A, all 
modes are completely frozen by the A S k term because, 
for all q , R\(q 2 ) is very large. Thus, Z k ^\ can be com¬ 
puted by the saddle-point method and it is then straight¬ 
forward to show that the presence of the last term in 
Eq. 0 leads to Ta [</>]— S'= <^]. On the contrary, when 
k = 0, the definition of R k implies that R k =o{q 2 ) = 0 
and the original model is recovered: Z k =o[J] = Z[J\ and 
r fe=0 [</] = F [</>], with T[4>] the usual Gibbs free energy or 
generating functional of one-particle-irreducible correla¬ 
tion functions. 

The exact RG equation for T*, reads [23l - l25| : 

d t T k [<j)\ = ^ f d t R k (\x - y\)G k , xy [4>], (4) 

^ J xy 

where t = ln(fc/A), and G ktXy [(f)\ is the field-dependent 
propagator: 

G ‘“(d 2> + ^)- 1 .rg»,W = j|A|^. (5) 

III. SCALE INVARIANCE 

We now discuss how scale invariance emerges in the 
NPRG formalism. We first consider a scale-invariant 


model described by an action S sca i- As we discuss in 
the following, the existence of such a model is not neces¬ 
sary for our proof but it is convenient to imagine that it 
exists to motivate the form of the expected Ward identity 
(WI) of scale invariance. If such a model existed, this WI 
would be obtained by performing the following infinitesi¬ 
mal change of variables ip(x) —>• (p(x) + e(D r f, + x l _ l d^)ip(x) 
in the functional integral, with D$ the scaling dimension 
of the field, usually written in terms of the anomalous di¬ 
mension r] as Dcf, = (d — 2 + rj)/2. Actually, the analysis 
of this model and of its WI faces both UV and infrared 
(IR) problems. In the IR regime, non-analyticities are 
present and, in the UV, it is difficult to control math¬ 
ematically the continuum limit: A —>• oo. Let us first 
discuss the IR aspect. Since ASfc[0] acts as an IR regu¬ 
lator, the Wilson RG offers a solution to the IR problem: 
the regularized model, Eq. CD, is not scale invariant even 
if the original model associated with S SC a\ were and, thus, 
all rg 0 ({pi}) are regular contrarily to r^ 0 ({pj}). The 
price to pay for regularity is the breaking of scale invari¬ 
ance that manifests itself through a modification of the 
WI [see [H, [13 for situations where R k breaks symme¬ 
tries]. By enlarging the space of cutoff functions R k to 
arbitrary functions R k (x,y) that are neither constrained 
to satisfy CD nor to be invariant under rotations and 
translations [as also done in m , this modified WI for 
scale invariance, obtained from Eqs. urn reads 


f (D* + D v + D R )R k (x,y) 

Jxy 

+ [ {D x + D^ix) 


sr k 


5R k (x,y) 

sr k 

5(f>(x) 


= 0 


( 6 ) 


where D x = x^d x and Dr = 2 d — 2is the scaling 
dimension of R k which implies that the field renormal¬ 
ization in Eq. CD behaves as Z k oc k T By constraining 
Rk(q) to be of the form JT]), Eq. (0 can be conveniently 
rewritten (following IH): 


d t r k [<t>\ = - f(D x + D+Mx) (7) 

Jx °<p{ x ) 


Introducing dimensionless and renormalized quantities 
(denoted with a tilde) 

x = k~ l x (8) 

(j>{x) = k D ' l ’(f>(x) (9) 

f fc [0] = r k [<t>\ (io) 

Eq. (JJD rewrites: 

<9tffc[</>] = 0. (11) 


Eq. (1111) means that an hypothetical scale-invariant ac¬ 
tion S/cai would lead, in its regularized and not scale- 
invariant version S^ai + A S k , to a RG flow where F[</>] 
would be at a fixed point r*[</] for all values of t: 

d t r*[4>] = o. 
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The very structure of the Wilson RG (or NPRG) also 
solves the UV problem. Actual models have a natural 
UV cutoff A (e.g. a lattice spacing) at which is defined 
their microscopic action S. The momentum integrals are 
therefore cut-off at A and are thus UV finite. At scales of 
order A (i.e., if we consider correlation functions, in the 
regime where at least one external momentum or k or the 
magnetic field in appropriate units, is of the order of A), 
the model is not scale invariant. In fact scale invariance is 
an emergent property that appears in the IR when some 
parameter (such as the temperature) has been fine-tuned. 
We call S'crit the corresponding action. In the RG for¬ 
malism, scale invariance emerges in the IR regime when 
integrating the RG flow starting at Ta[^>] = Scrit[^ = cj>] 
because Tk gets close to a fixed point for large negative t, 
that is, k A. As discussed above, the fixed-point con¬ 
dition coincides with the WI for scale invariance in pres¬ 
ence of a regulator, see Eqs. mnD. As a consequence, if 
the RG flow is attracted towards an IR fixed-point, then 
scale invariance emerges in the universal, long-distance 
regime 2 . We stress that this discussion does not rely 
on the actual existence of a well-defined continuum limit 
associated with a scale-invariant action S sca \ which is, 
per se, an interesting issue, but that we do not need to 
address in the present work. 

When the microscopic action is slightly different from 
the critical one (e.g. choosing a temperature slightly 
away from the critical one), the RG trajectory approaches 
the fixed point and then stays close to it for a long RG 
“time”, before departing. In this situation, the corre¬ 
lation length £ is finite but large and the WI is almost 
fulfilled for momenta A p £ _1 . This defines the 
critical regime of the theory. The closer the microscopic 
action is tuned to the critical one, the larger the correla¬ 
tion length £, and the better the WI is fulfilled. 

Let us now make two comments. First, when 
k 0 , d t Tk[4>] becomes negligible compared to any k- 
independent finite scale and Tk —> T. In this limit, 
Eqs. GED become the usual WI of scale invariance, as 
expected. Second, the above analysis shows that among 
the continuous infinity of solutions of the fixed point 
equation dtTk[4>] = 0 , only those that are regular for all 
fields are acceptable since they must be the limit when 
k — > 0 of the smooth evolution of Tfc[</>] from SAit [</>]• 
There is generically only a discrete, often finite, number 
of such physical fixed-points . 3 

A characteristic feature of the physical fixed-points 
is that the linearized flow around them has a discrete 
spectrum of eigenvalues from which some critical ex¬ 
ponents can be straightforwardly obtained jl5|. The 


2 A running anomalous dimension can be defined by r]k = 
—dt log Zk- It is only around the fixed point that //;, approaches 
a fixed-point value which is nothing else than //. 

3 Two well-known exceptions to this rule are the line of fixed points 
of the 0(2) model in d = 2 DT and the discrete infinity of 
(multicritical) ^-invariant fixed points in d = 2 Dd. T-jj . 


discrete character of the eigenperturbations around a 
fixed point has been studied intensively by perturba¬ 
tive means. As for the Wilson RG, it has been stud¬ 
ied in detail in the particular case of the O (IV) models 
in I34l - f36j . Although obtained within the derivative ex¬ 
pansion of the exact RG flow ©, its discrete character 
certainly remains true beyond this approximation. The 
eigenvalues A are obtained from the flow by substitut¬ 
ing Tfc[^] —> r*[^] + e exp(Af) 7 [ 0 ] and retaining only the 
O(e) terms. (With our definition of t, a relevant operator 
has a negative eigenvalue.) This leads to the eigenvalue 
problem 

A 7 [4>} = ! (D* + D+Mx) 3- 

Jx 0 <p{x) 

~\J [(£* + Dr)t(x - y)] Gt. 7W GW. (12) 

where 27 = {T, y , z, w}, G*[<j)} is the dimensionless renor¬ 
malized propagator at the fixed point: G* = (T*( 2 )+r ) -1 
and r(x) is the dimensionless inverse Fourier transform 
of r(q 2 /k 2 ) defined in Eq. ©. 

We conclude from the above discussion that regularity 
selects among all the fixed-point functionals T*[ 0 ] those 
that are physical, that is, that can be reached by an RG 
flow from a physical action ,S’ cr it and that have a discrete 
spectrum of eigenperturbations. 


IV. SPECIAL CONFORMAL 
TRANSFORMATIONS 


Let us now study conformal invariance by following 
the same method as above. As in the case of scale in¬ 
variance, to motivate the form of the WI we imagine 
that a conformally-invariant model, associated with an 
action jS con f[<p] written in terms of a primary field ip, 
exists. We put aside for now the problem of the ex¬ 
istence of this model since, as discussed below, we do 
not need that it actually exists for our proof. If such 
a model existed, the modified WI would follow from 
performing the infinitesimal change of variables p(x) —> 
ip(x) + e tl (x 2 d li - 2x^x v d v + 2 ax^)(p(x) in Eq. ©. By 
considering general cutoff functions as in Eq. © , we find 
that it reads: 

fjKZ- D x x, + Kl - D R y, 

+ £ {K ;-2D txMx) ^ pO, (13) 


with KJj = x 2 d* — 2x l _ l x v d*. By specializing Rk to 
functions of the form Eq. © and requiring again that 
Zk oc k~ v : Eq. (fl3l) can be rewritten 


0 = Kl4 , ]s Jy i *-2D, xMl )-A±- ) 


1 

2 


d t Rk(\x 


y I) ( x m + y m) Gk,xy 


(14) 
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Again, this identity boils down to the usual WI of con¬ 
formal invariance in the limit k —► 0 where R k —> 0. 

At any fixed point, the scaling dimension of ££[</>] is 
fixed by Eq. (fldll to be —1. We thus define £(([</>] = 
/c£^ [(/)]. Its flow equation reads: 


r xyu- 

d t t^\ - E£[0] = (D x + Djfc) 

Jx 0(p{x) 

- i J_[(D S + D R )r(x - y)\ G fc , a2 V “ ^ G k , m (15) 


where Yi^ 2 \<j>] is the second functional derivative of ££. 

An important property of £{( is that, at the fixed point, 
it is the integral of a density with no explicit dependence 
on x M . Indeed, observe that the left-hand-sides of Eqs. © 
and Cl can be interpreted as the action of the genera¬ 
tors of dilatations T> and conformal transformations 1C M 
on r fc . 


VT k 



{D x + D+Mx) 


ST k 

5<j){x) 


+ [ (D x + D v + D R )R k (x, y) 

J xy 


sr k 


SR k (x,y) 


(16) 


+ L {K ~ DHX “ + K ’- DxV ‘‘ }Rt{x ' v) Sl^-y)’ 

(17) 


that is, £ M = JC^T k . Similar expressions can be obtained 
for the generators of translations V M and rotations J lt y. 


VyT k = / d^x) 


sr k 

6(/)(x) 


0 % + dl)R k {x,y ) 


6T k 


SR k (x,y) 


(18) 


J^iXk = 


A<j)(x) 


ST k 

8(j){x) 


sr k 


{Xf,d x + yil dy)R k {x,y) 

, xy SR k (x,y) J 

-[p*+v\ 


(19) 


It can easily be checked that the generators satisfy the 
algebra of the conformal group. In particular, applying 
[Pn,K. v \ = 2(5^2? + 2 to a translation, rotation and 
dilatation invariant IA yields 


W = 0. (20) 

Thus, at the fixed point, £(!' is the integral of a density 
that does not have an explicit dependence on x. This 
density only depends on the field and its derivatives. This 
proof generalizes trivially to other scalar models. 


V. A SUFFICIENT CONDITION FOR 
CONFORMAL INVARIANCE 


Let us now consider a physical model at criticality. At 
the scale A, Ta = S = SAit and the model is neither con¬ 
formally nor scale invariant. However, when k <C A, the 
regularized model gets close to a fixed point and thus 
G k ~ G* and ££[^>] — [</>]. The key point of our 

proof is that at the fixed point, Eq. m is identical to 
m with 7 \(j)\ —» and A = — 1 although these two 

equations have different physical meanings. A sufficient 
condition for proving conformal invariance is therefore to 
show that there is no integrated vector eigenperturbation 
Vp = / V M of T*[4>] of scaling dimension Dy = —1. If no 
such eigenperturbation exists, the conformal WI is satis¬ 
fied, which means that the system is conformally invari¬ 
ant at criticality in the long distance regime. Moreover, 
the form of the conformal WI (ThEl) fixes the transforma¬ 
tion law fulfilled by <^>, which is the one of a primary 
field. 

To understand how conformal invariance is related 
with the scaling dimension of the vector eigenperturba- 
tions, it proves useful to consider three simple examples. 
The first one is the Ising model in d = 4. The fixed point 
being gaussian, the eigenvalues are trivially given by the 
canonical dimensions. By using the fact that £(( is Z 2 — 
and translation-invariant, see Eqs. (fl4l) and (l20l) . we find 
by inspection that the vector operator with lowest di¬ 
mension reads J <f> d^^dtp) 2 . It has therefore dimension 
+3. Note that there exists local vectors operators with 
lower scaling dimensions [(pd^cp, (p 3 d^cp and <pdfj,(p(d v (p) 2 )\. 
However these are total derivatives and are not associ¬ 
ated with integrated vector operators. In the absence of 
vector operator of dimension — 1 , we retrieve the well- 
known property that this model is conformally invariant 
at criticality in the long distance regime [H, [3. Using 
standard methods, we can compute the corrections to the 
scaling dimension of the vector operator /). (p d fi (p(d(p) 2 in 
a systematic expansion in e = 4 — d. We performed the 
calculation at one loop and found that the correction 
vanishes. The scaling dimension is therefore 3 + 0(e 2 ). 

This analysis can be extended to the 0(N) mod¬ 
els. In d = 4, there exists now two independent 
integrated vector operators, J (p a [d IJ ,(p a )(d v (p b ) 2 and 
j <p a {d v (p a ){d^<p b ){d v 4> b )^ with lowest scaling dimension 
3. As in the Ising case, there exist local operators with 
lower scaling dimensions, that are however total deriva¬ 
tives. Our sufficient condition is again fulfilled and we 
recover the well-known fact that these models are confor¬ 
mally invariant in the critical regime for d = 4. At one 
loop the degeneracy of the scaling dimensions is lifted 
and we obtain 3 + 0(e 2 ) and 3 — + 0(e 2 ). 

The third example involves a vector field A M (x) and is 
described by the (euclidean) action 


S = 



V 


duAfj ,) 2 + f (d.A ^) 2 


( 21 ) 


This model is interesting because it is scale invariant 
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FIG. 1: Possible behavior of the lowest eigenvalue Dy asso¬ 
ciated with a vector perturbation as a function of dimension. 
Left panel: (a) and (b) correspond to typical behavior, (c) 
to the exceptional case where Dv = — 1 right in d = 3. In 
the three cases, conformal invariance holds. Right panel: the 
shaded area represents a continuum of eigenvalues and the 
curve an eigenvalue Dv having a plateau at —1 around d = 3. 
In these cases, conformal invariance can be broken. 


but generically not conformally invariant, except for 
a = a c = (d — 4)/d [53, [HJ]. This situation can be under¬ 
stood in our context by considering the contraposition of 
our sufficient condition, which states that, assuming that 
is a primary field [38j, a necessary condition for hav¬ 
ing scale invariance without conformal invariance is that 
there exists an integrated vector operator with scaling di¬ 
mension — 1. It is actually easy to find such an operator 
(which is unique, up to a normalization): C f x A M (d u A u ). 
We can understand the particular case a = a c by an ex¬ 
plicit calculation of from Eq. (THT) . This shows 

that C = ad + 4 — d which, as expected, vanishes when 
a = a c . 

To conclude this section, we discuss the plausibility of 
not having conformal invariance in d = 3 for the O (TV) 
models at criticality. The only possibility would be to 
have a vector eigenperturbation with eigenvalue —1 in 
this dimension, see Fig. [I] This would mean that the 
d = 3 model has an integer critical exponent, a prop¬ 
erty that is highly improbable. Let us anyhow suppose 
that one of these eigenvalues crosses —1 right in d, = 3, 
as in curve (c) of Fig. [l] Then, for any dimension in¬ 
finitesimally smaller or larger than three, there would 
exist no eigenperturbation of dimension —1. The criti¬ 
cal system would exhibit conformal symmetry above and 
below d = 3. Since correlation functions of the criti¬ 
cal theory are expected to be continuous functions of d, 
we conclude that, even in this highly improbable situ¬ 
ation, the model would exhibit conformal invariance at 
criticality in d = 3. We are thus led to the more strin¬ 
gent necessary (but not sufficient) condition: for a critical 
model not to be conformally invariant, there must exist 
a vector perturbation of scaling dimension —1 in a finite 
interval of dimensions containing d = 3. This could hap¬ 
pen either because a discrete eigenvalue is independent of 
the dimension in some range of dimension around three 
or because there exists a continuum of eigenvalues, see 
Fig. m Such a behavior is, to say the least, not standard. 
To our knowledge, this has never been observed in any 
interacting model. 


The previous arguments are compelling but not mathe¬ 
matically rigorous. In particular, assuming that the the¬ 
ory can be properly defined in noninteger dimensions, 
which is standard but not obvious, it is hard to control 
the analytic structure of the critical exponents in d. In 
the next section, we give a proof in the physically impor¬ 
tant case of the Ising model in d = 3 which does not rely 
on such arguments. 


VI. PROOF IN THE ISING UNIVERSALITY 
CLASS 

We concentrate here on the Ising universality class and 
show that, in this case, the smallest eigenvalue Dy as¬ 
sociated with an integrated vector perturbation is larger 
than —1 for d < 4. Using our necessary condition, this 
proves that, for the Ising universality class, the critical 
regime is conformally invariant. 

For simplicity, we consider a lattice version of the 
Ginzburg-Landau model whose dynamics is described by 
the hamiltonian (or action) 

S =+^2U(<p(i)). (22) 

<U> * 

where the index i labels the lattice sites, the ip{i) take 
values in the real domain and U{ip) is an even function 
that diverges for \ip\ —> oo. We choose here a hypercubic 
lattice with lattice spacing a. The original Ising model 
can be recovered by considering a potential U {ip) strongly 
peaked around ip = ±1 but the Ginzburg-Landau model 
which is in the Ising universality class for a generic po¬ 
tential is more convenient for what follows. 

The use of Ginzburg-Landau model has another ad¬ 
vantage. In the case of a quartic potential, 

£%(*)) = yP 2 (*) + ( 23 ) 

and for d < 4, the model is super-renormalizable. Its ul¬ 
traviolet behaviour is therefore controlled by a Gaussian 
fixed-point. In this case, the existence of a controlled 
scaling limit seems to be under control even if, to the 
best of our knowledge, there is no mathematical proof 
of its existence [39j. To compare different models in the 
Ising universality class we assume below that this scal¬ 
ing limit does exist in the following precise sense. Let us 
consider two local operators Oi{i, a) and €> 2 ( 1 , a). Let us 
also introduce (i) a smooth interpolating field c/) lnteip {x) 
with x GM. d that coincides on the lattice points with cf>{i) 
and (ii) two interpolating operators {x, a) defined 

by: 0) nterp = Oi{(f>' nterp ). Of course the construction 
of these interpolating operators is not unique. We now 
consider the particular case where 0\ and O 2 are such 
that 0) nterp —>. ( 9 “ tcl P w hen a —I 0. (We notice that 
this limit is independent of the choice of interpolation 
used to define <(> mterp ). Our assumption, that we call for 
short “scaling limit”, is that there exists a multiplicative 
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factor Zo(a) depending on the lattice spacing such that 
the correlation functions of the operators 0\{x,a) and 
Zo(a)0 2 {x,a) are the same for distances much larger 
than a: 

(0 1 (x,a)0 3 (y 3 )...0 n (y n )) 

~ Z 0 (a)(0 2 (x,a)0 3 (y 3 )... O n (y n )) (24) 


where d^ is a lattice discretization of the partial deriva¬ 
tive. (Note that the first three operators are total deriva¬ 
tives.) In dimension d < 4, although we do not know the 
explicit form of the eigenoperators, we can, in principle, 
decompose any vector operator on the basis and 

generically, there is a nonvanishing overlap with each of 
the eigenoperators: 


where {0 3 ,... ,O n } are arbitrary local operators and 
where the equivalence occurs for a -C min {|?/3 — 
x\,... ,\y n — x\}. This hypothesis is a pre-requisite of 
all Monte Carlo simulations, and is, of course, satisfied 
to all orders of perturbation theory in any renormalizable 
theory. We assume here that it is also valid nonpertur- 
batively. 

Our strategy is to study correlation functions of local 
vector operators V^(x) and use their critical behavior to 
find a bound on the scaling dimension of the integrated 
operator We rightaway mention two 

difficulties. 

First, there arc local vector operators that are total 
derivatives and which are therefore not associated with 
an integrated one. As discussed before, the vector oper¬ 
ator d^p 2 ) is such an operator. Note that its scaling 
dimension near d = 4 is lower than that of the vec¬ 
tor operators which are not total derivatives (such as 

Second, operators with the same quantum numbers 
typically mix together in the calculation of correlation 
functions. As a consequence, the critical behavior of a 
two-point function of some vector operator is governed 
by the lowest scaling dimension of the class of opera¬ 
tors with which it mixes. To be more precise, let us call 
the local vector eigenoperator of scaling transfor¬ 
mations with scaling dimension D^ 1 (ordered such that 

Dy 1 < Dy' 1 < •••). The associated two-point correla¬ 
tion function behaves, in the critical regime, as: 

(V'-’Wt’M.-;-(25) 

\x-y\ 2U v 

where the subscript c indicates connected correlation 
functions, defined as: 

(X(x)Y(y)) c = (X(x)Y(y)} - (X(x))(Y(y)) (26) 


V/* = I>" v M n) - ( 31 ) 

n 

As a consequence, the critical regime of the two-point 
correlation function is dominated by the smallest critical 
dimension 

<V^(*)V M (i,)> c ~--(32) 

\x — y\ v 

We stress that the list of quantum numbers must include 
those associated with lattice isometries. For example, 
we require scalar (respectively vector) operators defined 
on the lattice to be even (respectively odd) under parity 
transformations. 

The proof is organized as follows. Using Lebowitz in¬ 
equalities Hc],[4lj], we derive a lower bound for from 
which follows a lower bound for the scaling dimension of 
the integrated vector operators. The proof that the scal¬ 
ing dimension Dy of any integrated vector operator is 
different from —1 for d < 4 is then a direct consequence 
of this bound. 

As a first step, we derive a bound for the correlation 
function ( t p 2 ( x)ip 2 (y)) r ,. We use here the Lebowitz in¬ 
equalities [40| which state that, considering two decou¬ 
pled copies of the ferromagnetic system (that we note tp 
and p r ), both described by the action (l22|) . and consid¬ 
ering two sets A and B of lattice points, 


( n [^(*) + - ¥>'l7)]> < 

i&A j£B 

( n [pW + v'WMj) + <p'ti )]> > 

i&A,jeB 

di^w+^wixn^+^o-)])- 

i£A j£B 


(33) 


(34) 


and where an appropriate normalization has been chosen. 
In d = 4, the eigenproblem can be solved and the scaling 
dimensions are the canonical dimensions. In particular, 
the eigenoperators with lowest scaling dimensions are: 


V£°’ d = 4) oc d ^ 2 
V< M=4) cx 
Vf d=4) oc d^d ^) 2 
Vf d=4) cx (d^ 2 )(d ^) 2 


7-)(0, d— 4) 
U V 

= 3 

(27) 

7-)(l»^=4) 

U V 

= 5 

(28) 

r-)(2,d=4) 

U V 

= 5 

(29) 

p.(3,d—4) 

U V 

= 7 

(30) 


In particular, this implies that: 

{(<p(x) + p’(x)) 2 {ip(y) - p\y)f) < 

((p(x) + p'(x)) 2 )((p(y) - p'{y)) 2 ). 

(35) 

Expanding the binomials, we readily obtain the following 
identity: 

{ip 2 (x)<p 2 (y)) c < 2G 2 (x - y) (36) 

where we have used the fact that the average of an odd 
number of fields vanishes for temperatures higher than 
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(or equal to) the critical temperature. This implies that 
the connected correlation function (<p 2 (x)tp 2 (y)) c cannot 
decrease more slowly than the square of the propagator 
at long distances. This inequality can be generalized to 
arbitrary even powers of the Helds: 

0 < (ip m (x)<p n (y)) c < C G 2 (x — y) (37) 

where C is a positive constant (that depends on m and 
n) as shown in the Appendix [A] 

In the critical regime, scale invariance implies that con¬ 
nected two-point correlation functions behave as power- 
laws. In particular: 

(y-Wrt)), ~ |x _ (38) 

with A a positive constant (see, for example 0 )- The 
inequality (1571) implies that > d — 2 + r]. We can then 
deduce the asymptotic behavior of the matrix of second 
derivatives of this correlation function: 

~ir—jraraK'- 

1 y w , ( 39 ) 

, ^ (z - wz - 2/M 

+ c {^y) 2 )■ 

where B and C are some constants. 

We now consider two local vector operators that are 
the product of one power of d^ip(x) and an odd (finite) 
number of fields evaluated at points in a finite neighbor¬ 
hood of x: 

1 m— 1 

W^ta:) = -[d^ix)] Y v( x + s e ^)- ( 40 ) 

S —±1 2=1 
.. n—1 

= -\d^{x)\ Y n v{x + sef ] ). (41) 

S —dzl 2 = 1 

where and e\■ 2 ' > are some constant lattice vectors 4 . 

The operators W^\x) and wj?\x) are, up to a multi¬ 
plicative constant, other discretizations of, respectively, 
the operator d*((p m (x)) and d^(ip n (x)). According to the 
assumption of the existence of the scaling limit, Eq. ([24]) . 
the connected correlation function (W^ (x)W^ (y)) c 
has the same asymptotic behavior as in Eq. (1391) up to 
a multiplicative factor depending on the lattice spacing. 
Indeed, when \x — y\ is much larger than the lattice spac¬ 
ing a, the vectors e- 1 ^ and e- 2) can be neglected in (1701) 
and (1711) and the local operators and W^ 2 '* are then 
proportional to and d u (<p n ) respectively as ex¬ 

plained before. 


4 Since , defined in Eq. m is odd under parity, it is important 
in what follows to consider only vector operators that are also 
odd. This is the reason why the sum over s is necessary in the 
definitions (Hot and S3- 


Now, any local vector operator on the lattice is a lin¬ 
ear combination of vector operators of the form (17U1) . For 
instance, a discretization of the operator d ll ((t> 2 )(d v (l>) 2 is 
given by: 

{4>{x + A) - <t>{x - A)) Y x + c, )~ <l>( x ~ v)f 

(42) 

where the sum runs over all the nearest neighbors of x. 

Using the triangular inequality, we conclude that: 

l<r,MK(v)>cl < \ x _f£ d _^ r (43) 

where is a positive constant. Using Eq. (l32ll this 
implies that, for all n: 

£>v n) > d - 1 + Tj (44) 

We conclude that the scaling dimension Dy = Dy — d of 
any integrated vector operator is not smaller than — 1 + 77. 5 
Using the unitarity of the Minkowskian </> 4 theory, one 
can prove that 77 > 0 [H] for the Ising universality class. 
Moreover, an interacting massless theory such as the crit¬ 
ical Ising model for d < 4, has a non-zero rj (43j |. As a 
consequence, our necessary condition is fulfilled and we 
conclude that scale invariance implies conformal invari¬ 
ance for the Ising universality class for all d < 4. 

VII. CONCLUSIONS 

Let us now point out some directions of research for 
the future. It is clear that the condition of conformal in¬ 
variance (fl4l) can be straightforwardly extended to other 
theories (involving scalar, fermionic or vector fields) and 
it would be interesting to conclude on the fate of confor¬ 
mal invariance in this wider class of models. In these sys¬ 
tems, it is much more difficult to find rigorous bounds on 
correlation functions (that would generalize the Lebowitz 
inequalities). It would then prove useful to approach the 
problem by computing the scaling behavior of vector op¬ 
erators by Monte-Carlo simulations. 

Another promising line of investigation consists in 
making use of the conformal invariance in the Wilson 
framework to perform actual calculations of universal 
quantities. On the one hand, and in the best case, this 
would lead to closed and numerically tractable equations 
for the critical exponents. On the other hand, the ap¬ 
proximation schemes currently used for solving the Wil¬ 
son RG flow equation being incompatible with exact con¬ 
formal invariance, we can expect that constraining them 
to be conformally invariant at the fixed point would im¬ 
prove their accuracy. 


5 This bound also applies for correlation functions of more general 
operators. 









Notice finally that, at first sight, our approach could 
seem similar to the one based on the energy-momentum 
tensor and on the analysis of the virial current. This is 
not the case although there is perhaps a relationship be¬ 
tween the two. E(( is a functional of <f> and not of it 
is built from T& and not from S'; what matters is that 
its density vanishes up to a surface term and not that 
it is conserved. Moreover, as we already explained, we 
only deal with a regularized theory which enables us to 
consider only the analytic candidates for E^ 1 * contrary 
to what should be done in a nonregularized theory. In 
any case, a clarification of the relation between the two 
approaches would be welcome. In this respect, our proof 
of the non-existence of local vector operator of scaling di¬ 
mension d— 1 (conserved or not) might be of interest also 
when applied to a hypothetical conserved virial current. 
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Appendix A: Bound for correlation functions 

(<p m {x)tp n {y)) 

In this appendix, we derive bounds on the correlation 
functions ((p m (x)ip n (y)) c , in the symmetric phase (T > 
T c ), with to and n arbitrary integers with the same parity 
(otherwise the correlation function vanishes). 

We want to show that: 

(ip a (x)tp b (y)) c < C\G(x - y) for odd a, b (Al) 

(tp a (x)<p b (y)) c < C 2 G 2 (x — y) for even a, b (A2) 

where C\ and C 2 are some strictly positive constants and 

G{x ~ y) = {<p(x)p(y)). (A3) 

Note that for odd a and b, the connected and discon¬ 
nected correlation functions are equal, see Eq. liTOl 

Property ED is obvious for a = b = 1. The proof 
of (IA2I) for {a = 2,b = 2} is presented in the core of 
the article. For general a and b, the proof is made by 
induction. Assuming that the inequalities (ED and (IA2I) 
are fulfilled for {a < to, b < n} \ {a = m, b = n} , we have 
to prove that these properties are also valid for a = m 
and b = n. 


We first consider the case where m and n are even. 
Using the Lebowitz inequality [4(| [see Eq. ®] 

(\p(x) + tp'(x)] m [ip{y) - p\y)] n ) < 

Mx) + y'{x)r)My)-v\y)r) 

(A4) 

as well as translation invariance and the binomial expan¬ 
sion, we obtain: 


ED- 1 )' 


a—0 b— 0 


m \ In 


(p a (x)p b (y)) c (p m - a (x)p n - b (y)) 


( l P a ( 0 )) (<p ( 0 )) {(p m ~ a (x)ip n ~ b (y)) c 

(p a ( X )p b (y)) c (p m - a m(p n - b m 


< 0 


(AS) 


Writing explicitly the terms with a £ {0, m} and b £ 
{ 0 ,n}, we get the following bound: 

2{<p m (x)ip n (y)) c < E E 

a£{l,3,--- ,m—1} 6£{ 1,3, — ,n—1} 

(p a (x)p b (y)) c (p m - a (x)p n ~ b ( y )) c 

(A6) 

where we have used the fact that connected correlation 
functions, as well as (</? m ( 0 )) are non-negative [44|,|45| to 
eliminate the terms with even b. We observe that there 
appears in the right-hand side only connected correlation 
functions with at most to — 1 powers of (p(x) and at most 
n — 1 powers of p{y). By hypothesis, properties (IA1I) and 
m are thus valid for these correlation functions. Fur¬ 
thermore, we observe that the right-hand side involves a 
sum of terms that are a product of two correlation func¬ 
tions with odd powers of the fields. In both cases, by 
using properties ED and (1A2I) , this quantity is bounded 
by some positive constant times G 2 (x — y). 6 This con¬ 
cludes the proof of the induction hypothesis in the case 
of even to and n. 

We now turn to the case where to and n are odd. We 
now make use of the Lebowitz inequality 

Mx) + p'{x)r- l [p{x) - p'(x)My) - p'(y)D < 

Mx) + ip'Wr-^Mx) - p'(x)My) - ip'{y)] n ) 

(A7) 



which is of interest if m > 1 (we can obviously derive 
a similar inequality with {m,x} {n,y} which can be 

applied in the case m = 1). We again use the binomial 


We use here the fact that, on the lattice, G is bounded by some 
positive constant ( G(r ) < D ). Consequently, we can further use 
the bound G 2n (r) < CG 2 (r) where n is a positive integer and C 
a positive constant. 
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expansion and the positivity of (connected and discon¬ 
nected) correlation functions to obtain the following in¬ 
equality: 


m — 1\ f n 


(<p m (x)<p n (y)) c < 

E E 

a6{0,2, -- ,m- 1} be{l,3,--- ,n- 1} 

(v o (0))(v m - a - 1 (0))^(*)v 6 (|/)) c ^ n - 6 (0)) (A8) 

+ E E 

aG{l,3,-" ,m-2} be{0,2,---,ra} v ' 

( t p a+1 { x ) ( p b (y))( i p m ~ 1 ~ a ( x ) i p n ~ b (y)) 


tion functions, one with even and one with odd powers of 
the fields, or the product of a correlation function with 
odd powers of the fields and a positive constant. In all 
cases, the correlation functions that appear in the right- 
hand side fulfill the conditions of our hypothesis. We 
therefore conclude that (ip m (x)ip n (y)) c satisfies property 
ED for to and n odd (see Footnote 6). This concludes 
the proof of the induction hypothesis. 


Using the fact that the property ED is valid for a = 
b = 1, it is easy to check, by applying several times the 
induction property, that ED and (IA2I) are valid for any 
a and b. 


The first term involves the product of a correlation func¬ 
tion with odd powers of the fields and a positive constant. 
The second sum involves either a product of two correla- 
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